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Sedimentary  organic  matter  from  the  Here’s  Your  Chance  (HYC)  Pb-Zn-Ag  deposit  (McArthur  Basin, 
Northern  Territory,  Australia)  displays  increased  thermal  maturity  compared  to  nearby  non-mineralised 
sediments.  Micro-scale  sealed  vessel  pyrolysis  (MSSVpy)  of  an  immature,  organic  rich  sediment  from  the 
host  Barney  Creek  Formation  (BCF)  was  used  to  simulate  the  thermal  maturation  of  OM  from  the  HYC 
deposit,  and  to  assess  the  effect  of  sulfide  minerals  on  organic  maturation  processes.  MSSVpy  at  increas¬ 
ing  temperatures  (300,  330  and  360  °C)  resulted  in  increased  methylphenanthrene  maturity  ratios  which 
were  within  the  range  reported  for  bitumen  extracted  from  HYC  sediments.  The  methylphenanthrene 
index  ratio  from  MSSVpy  of  the  BCF  sample  was  lower  than  in  HYC,  due  to  a  reduced  proportion  of  meth¬ 
ylated  phenanthrenes.  Polycyclic  aromatic  hydrocarbons  with  four  to  six  rings  were  produced  in  increas¬ 
ing  abundance  as  pyrolysis  temperature  increased,  although  they  did  not  approach  the  levels  reported 
from  HYC.  Addition  of  lead  sulfide,  zinc  sulfide  and  pyrite  to  the  MSSVpy  experiments  resulted  in  a 
reduced  response  of  the  methylphenanthrene  maturity  parameters,  possibly  due  to  retardation  of 
methyl-shift  and  transmethylation  reactions. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Organic  matter  (OM)  is  thought  to  play  an  important  role  in  the 
formation  of  many  commercially  viable  ore  bodies,  hence  there  is 
significant  interest  in  using  organic  geochemical  techniques  to 
investigate  the  formation  history  of  ore  deposits  (Greenwood 
et  al.,  2013).  OM  can  be  significantly  altered  by  interaction  with 
ore-forming  fluids.  For  example,  major  hydrothermal  events  asso¬ 
ciated  with  the  formation  of  sedimentary-exhalative  (SEDEX)  ore 
deposits  have  caused  a  marked  increase  in  the  thermal  maturity 
of  associated  sedimentary  OM  (Gize,  1999;  Logan  et  al.,  2001; 
Williford  et  al.,  2011).  Oxidised  fluids  which  formed  the  Kupfers- 
chiefer  deposit  contributed  to  increased  pristane/phytane  ratios 
(Large  and  Gize,  1996)  and  lower  abundances  of  methylated 
aromatics  (Piittmann  et  al.,  1989).  The  high  abundance  of  sulfide 
minerals  in  such  environments  may  be  a  major  influence  on  the 
behaviour  of  OM  during  alteration,  such  that  traditional 
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assumptions  and  interpretations  largely  based  on  petroleum 
organic  geochemistry  may  not  apply  (Gize,  1999).  Although  sulfide 
minerals  have  been  found  to  catalyse  abiotic  organic  reactions 
important  to  the  early  origins  of  life  (Cody  et  al.,  2004;  Cleaves 
et  al.,  2012)  their  effect  on  the  thermal  maturation  of  OM  has  been 
largely  unexplored. 

A  wide  range  of  molecular  parameters  have  been  developed  to 
assess  the  thermal  maturity  of  source  rocks  and  petroleum  hydro¬ 
carbons,  summarised  by  Peters  et  al.  (2005).  Parameters  based  on 
polycyclic  aromatic  hydrocarbons  (PAHs)  such  as  naphthalene 
and  phenanthrene  are  widely  used  due  to  the  ubiquity  of  these 
compounds  in  sediments  and  their  high  thermal  stability 
compared  to  aliphatic  biomarkers  (Price,  1993).  Methylated  PAHs 
readily  undergo  isomerisation  reactions  during  maturation,  with 
the  (S-substituted  isomers  being  more  thermodynamically  stable 
than  the  ot-isomers  (Nabbefeld  et  al.,  2010;  Szczerba  and 
Rospondek,  2010;  Dick  et  al.,  2013).  A  number  of  parameters  have 
been  developed  based  on  the  relative  abundance  of  methylphe¬ 
nanthrene  (MP)  isomers,  including  the  methylphenanthrene  ratio 
MPR=2-MP/1-M  (Radke  et  al  1982b)  and  the  methylphenanth¬ 
rene  indices  MPi-1  =  1.5  (2-MP  +  3-MP)/(P  +  9-MP  +  1-MP)  (Radke 
et  al.,  1982a)  and  MPI-3  =  (2-MP  +  3-MP)/(9-MP  +  1-MP) 
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(Garrigues  et  al.,  1988).  Catalysis  by  clay  mineral  surfaces  has  been 
shown  in  different  circumstances  to  promote  methylation  or 
demethylation  of  PAHs  in  shales  and  crude  oils  (Alexander  et  al., 
1995;  van  Aarssen  et  al.,  1999)  and  thus  can  influence  the  values 
of  aromatic  maturity  parameters  (Nabbefeld  et  al.,  2010).  Larger 
PAHs  with  four  or  more  aromatic  rings  are  not  generally  used  in 
maturity  parameters,  although  they  may  be  produced  by  burial 
metamorphism  (Brocks  et  al.,  2003)  and  have  been  detected  in 
high  abundances  in  submarine  hydrothermal  petroleum  (Kawka 
and  Simoneit,  1990;  Simoneit,  1994;  Simoneit  et  al.,  2004).  PAHs 
have  been  detected  in  a  range  of  ore  deposits  including  sedi¬ 
ment-hosted  lead-zinc  (Chen  et  al.,  2003;  Williford  et  al.,  2011), 
mercury  (Blumer,  1975;  Wise  et  al.,  1986),  copper-silver 
(Puttmann  et  al.,  1989;  Rospondek  et  al.,  2007,  2009)  and  strata- 
bound  copper  sulfide  (Rieger  et  al.,  2008). 

Micro-scale  sealed  vessel  pyrolysis  (MSSVpy)  is  a  closed  system 
pyrolysis  technique  in  which  relatively  small  (mg)  amounts  of 
organic  material  are  heated  in  a  sealed  glass  tube.  MSSVpy  has 
been  widely  applied  to  simulate  the  catagenetic  reactions  of 
hydrocarbons  and  study  the  kinetics  of  petroleum  generation 
(e.g.  Horsfield  et  al.,  1989;  Schenk  and  Horsfield,  1993; 
Ladjavardi  et  al.,  2013),  as  well  as  to  characterise  naturally  occur¬ 
ring  OM  (Greenwood  et  al.,  2006;  Berwick  et  al.,  2010;  Lockhart 
et  al.,  2013)  and  investigate  the  formation  of  specific  organic  prod¬ 
ucts  during  diagenesis  (Sinninghe  Damste  et  al.,  1998).  Another 
common  technique  for  the  artificial  maturation  of  OM  is  hydrous 
pyrolysis,  in  which  crushed  or  powdered  rock  is  heated  with  water 
in  sealed  vessels  (Lewan  et  al.,  1979).  Compared  to  MSSVpy, 
hydrous  pyrolysis  more  holistically  replicates  the  conditions  of 
natural  organic  maturation  due  to  the  presence  of  water  and  the 
control  of  pressure,  fluid  composition  and  redox  state  (e.g. 
Landais  et  al.,  1994;  Lewan,  1997;  Seewald,  2001).  However  the 
high  pressures  generated  during  hydrous  pyrolysis  require  specia¬ 
lised  laboratory  equipment,  while  MSSVpy  is  easier  to  employ  and 
generates  similar  products  to  natural  organic  maturation 
(Horsfield  et  al.,  1989).  Here  we  employ  MSSVpy  to  simulate  the 
thermal  maturation  of  OM  from  the  Here’s  Your  Chance  (HYC) 
Pb-Zn-Ag  deposit  and  to  assess  the  effect  of  high  abundances  of 
sulfide  minerals  on  aromatic  maturity  parameters. 


2.  Materials  and  methods 

2.1.  Sample  site  and  description 

The  Barney  Creek  Formation  (BCF)  is  a  1640  ±3  Ma  (Page  and 
Sweet,  1998)  marine  shale  deposited  in  the  McArthur  Basin 
(Northern  Territory,  Australia)  under  reduced,  deep  water  condi¬ 
tions  (Bull,  1998).  Detailed  descriptions  are  presented  by  Brocks 
and  Schaeffer  (2008)  and  Greenwood  et  al.  (2013).  The  BCF  is  host 
to  the  Here’s  Your  Chance  (HYC)  Pb-Zn-Ag  deposit,  the  largest  of 
the  Proterozoic  deposits  of  the  North  Australian  craton  which  is 
the  world’s  richest  Pb-Zn  mineral  province  (Huston  et  al.,  2006). 
The  base  metals  in  HYC  are  widely  believed  to  have  been  deposited 
by  a  hydrothermal  brine  with  a  temperature  of  100-200  °C 
(Huston  et  al.,  2006;  Williford  et  al.,  2011),  in  a  tectonically  con¬ 
trolled  sub-basin  (McGoldrick  et  al.,  2010).  A  number  of  studies 
have  investigated  the  organic  geochemistry  of  HYC  (Logan  et  al., 
2001;  Chen  et  al.,  2003;  Williford  et  al.,  2011)  and  have  found  evi¬ 
dence  of  significant  hydrothermal  alteration  of  OM  in  mineralised 
sections.  Mineralised  sections  of  the  deposit  contain  abundant 
fine-grained  pyrite,  sphalerite  and  galena  (e.g.  Large  et  al.,  1998; 
Ireland  et  al.,  2004). 

The  sample  used  in  this  study  (referred  to  as  the  BCF  sample) 
was  taken  from  the  Clyde  Sub-basin,  from  drill  core  GR-7  at  depth 
94.69-94.81  m.  The  Glyde  Sub-basin  is  located  80  km  south  of  HYC 


and  is  stratigraphically  equivalent  and  morphologically  similar  to 
the  deposit  (Davidson  and  Dashlooty,  1993;  McGoldrick  et  al„ 

2010).  Sediments  from  the  Glyde  region  often  have  low  thermal 
maturity  for  samples  of  Paleoproterozoic  age,  with  equivalent  vitr- 
inite  reflectance  of  0.6-0.8%  calculated  from  MPI-1  (Crick,  1992), 
and  contain  abundant  biomarkers  of  phototrophic  sulfur  bacteria 
including  intact  carotenoid  pigments  (Brocks  et  al.,  2005;  Brocks 
and  Schaeffer,  2008).  Organic  matter  preserved  within  the  kero- 
gen-mineral  matrix  of  mineralised  HYC  samples  has  shown  evi¬ 
dence  for  the  presence  of  these  bacteria  during  ore  formation 
(Holman  et  al.,  2014).  The  BCF  sample  was  deposited  in  a  similar 
environment  to  HYC  and  displays  relatively  low  thermal  maturity. 
While  the  sample  is  not  an  exact  match  of  the  HYC  sediments  prior 
to  mineralisation,  it  is  a  suitable  analogue  to  evaluate  the  effect  of 
sulfide  minerals  on  the  organic  maturation  processes  that  occurred 
during  the  formation  of  the  deposit. 

2.2.  Sample  preparation 

The  BCF  sample  had  been  stored  long  term  in  a  Whirl-Pak  plas¬ 
tic  sample  bag.  Such  plastic  bags  are  potential  sources  of  hydrocar¬ 
bon  contamination  including  branched  alkanes  with  quaternary 
carbons  (BAQCs;  e.g.  Grosjean  and  Logan,  2007;  Greenwood 
et  al.,  2009).  Before  the  sample  was  powdered,  rock  pieces  were 
ultrasonically  extracted  (30  min)  with  a  small  amount  of  Milli-Q 
purified  water  to  ablate  the  surfaces  of  the  sample.  The  water 
and  fine  suspended  particles  were  twice  shaken  with  DCM  in  a  sep¬ 
arating  funnel.  The  organic  layer  was  reduced  with  a  rotary  evap¬ 
orator  and  passed  through  a  Pasteur  pipette  packed  with  MgS04 
(Fisher  Scientific,  pre-rinsed  with  DCM).  Three  small  pieces 
(lxl  cm)  were  cut  from  the  sample  bag  with  scissors  (rinsed  with 
DCM)  and  ultrasonically  extracted  for  15  min  in  DCM.  Both 
extracts  were  analysed  by  GC-MS  (instrumental  conditions 
detailed  in  Holman  et  al.,  2012).  Analysis  of  the  m/z  99-85 
extracted  ion  chromatogram  from  the  sample  bag  extract  revealed 
a  series  of  3,3-diethylalkanes  with  21-35  carbon  atoms,  which  are 
known  contaminants  from  Whirl-Pak  sample  bags  (Grosjean  and 
Logan,  2007).  These  compounds  were  not  detected  in  the  surface 
washings  from  the  BCF  sample. 

The  sample  was  ground  to  a  fine  powder  in  a  Rock  Labs  ring  mill 
with  a  zirconia  head.  The  mill  was  cleaned  before  use  by  repeated 
grinding  of  acid  washed,  annealed  sand.  The  final  batch  of  ground 
sand  was  extracted  with  DCM  and  GC-MS  analysis  showed  no  sig¬ 
nificant  contamination. 

2.3.  Pure  metal  sulfides 

Zinc  sulfide  (Purum,  Riedel-de  Haen,  Sigma  Aldrich)  was 
received  as  a  fine  white  powder.  Lead  sulfide  (99.9%,  Aldrich) 
was  in  the  form  of  aggregates  that  were  crushed  to  a  fine  powder 
with  a  mortar  and  pestle.  Pyrite  (Alfa  Aesar)  was  received  as  grains 
approximately  1.5-5  mm  in  length,  which  were  powdered  in  a 
Rock-Labs  ring  mill  with  a  zirconia  head.  Sulfides  were  cleaned 
by  ultrasonic  extraction  (x3)  with  DCM.  Analysis  of  the  final 
extract  by  GC-MS  showed  no  significant  contaminants. 

To  examine  the  size  distribution  of  the  sulfides,  particles  were 
suspended  in  deionized  water,  transferred  to  microscope  slides 
and  imaged  on  an  optical  microscope  under  transmitted  light. 
The  sulfides  consisted  of  fine  particles,  often  sub-micrometre  in 
size,  commonly  aggregated  into  clusters  ranging  from  ten  to  sev¬ 
eral  hundred  micrometres.  Previous  investigations  of  sulfide  min¬ 
eral  textures  in  the  HYC  deposit  have  found  that  sulfide  minerals 
occur  in  grains  of  1-10  pm,  aggregated  into  clusters  of  20- 
200  pm  (Eldridge  et  al.,  1993;  Ireland  et  al.,  2004).  The  pure  sul¬ 
fides  used  in  this  study  are  therefore  an  adequate  approximation 
of  the  sulfide  minerals  occurring  in  HYC. 
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Table  1 

Precise  weights  of  BCF  sample  and  metal  sulfide  for  the  mixtures  sub-sampled  for 
MSSVpy  experiments. 


Sample  name 

BCF  sample  (mg) 

Metal  sulfide  (mg) 

Sulfide  (wt%) 

5%  PbS 

961.45 

54.04 

5.32 

20%  PbS 

865.98 

216.57 

20.01 

5%  ZnS 

968.86 

50.96 

5.00 

20%  ZnS 

783.12 

196.76 

20.08 

5%  Py 

954.88 

49.26 

4.91 

20%  Py 

796.24 

196.90 

19.83 

Mixtures  of  BCF  sample  with  approximately  5  wt%  and  20  wt% 
metal  sulfides  were  prepared  and  homogenised  with  a  mortar 
and  pestle.  Precise  weights  for  all  mixtures  are  listed  in  Table  1. 

2.4.  Micro-scale  sealed  vessel  pyrolysis  and  GC-MS  analysis 

MSSVpy  followed  a  previously  described  procedure 
(Greenwood  et  al.,  2006;  Berwick  et  ah,  2010).  Extractable  organic 
matter  was  not  removed  from  the  sample  prior  to  pyrolysis.  Glass 
tubes  (5  cm  x  5  mm  i.d.)  and  glass  beads  (Alltech  60/80)  were 
cleaned  by  heating  overnight  at  300  °C  before  use.  Typically  6- 
10  mg  of  sample  was  loaded  into  the  centre  of  the  tube  and  the 
void  space  above  and  below  the  sample  was  filled  with  glass  beads. 
The  tubes  were  flame  sealed  before  heating  in  an  oven  for  72  h  at 
either  300  °C,  330  °C  or  360  °C.  The  tubes  were  then  placed  in  a  ser¬ 
ies  II  MSSV  injector  (held  at  300  °C)  mounted  on  a  GC-MS  (Agilent 
6890  GC  coupled  to  an  Agilent  5973  MSD)  and  cracked  with  a  plun¬ 
ger.  Volatile  products  were  transferred  to  the  GC  column  (Agilent 
J&W  DB5-MS,  60  m  x  0.25  mm  i.d.,  0.25  pm  film  thickness)  with 
He  carrier  gas  at  a  constant  pressure  of  9  psi,  and  a  split  ratio  of 
20:1.  The  oven  was  held  at  40  °C  for  2  min  while  the  column  was 
immersed  in  liquid  nitrogen  to  cryogenically  focus  the  products. 
The  oven  temperature  was  increased  from  40  °C  to  320  °C  at 
4  °C/min  then  held  at  320  °C  for  35  min.  Mass  spectra  were 


acquired  in  full  scan  mode  (m/z  25-550)  with  an  ionisation  energy 
of  70  eV.  All  samples  were  prepared  and  analysed  in  duplicate. 
Product  assignments  were  based  on  GC  retention  time  and  mass 
spectral  comparison  with  standards  or  library  data. 

2.5.  Correction  factors 

Due  to  the  online  nature  of  MSSVpy  the  pyrolysates  cannot  be 
separated  into  saturate  and  aromatic  fractions  prior  to  GC-MS 
analysis.  The  GC  chromatograms  produced  in  this  study  contain 
products  which  fully  or  partially  co-elute  and  accurate  integration 
of  the  total  ion  chromatogram  (TIC)  was  often  difficult  for  com¬ 
pounds  of  low  abundance  (see  Fig.  1  for  a  representative  chromato¬ 
gram).  For  this  reason  aromatic  compounds  were  integrated  from 
extracted  ion  chromatograms  (EIC)  of  diagnostic  ions.  The  selected 
ion  makes  up  a  different  proportion  of  the  total  mass  spectrum  for 
each  compound;  hence  ratios  calculated  from  EIC  cannot  be 
directly  compared  to  those  from  TIC.  Correction  factors  were  thus 
determined  for  each  aromatic  compound  of  interest  by  calculating 
the  ratio  of  peak  area  from  EIC  to  TIC  from  previous  samples  run  on 
the  GC-MS  instrument.  Only  compounds  free  from  co-elution  were 
included  in  the  calculations,  and  a  range  of  peak  areas  were  chosen 
to  ensure  linearity  of  the  calibration.  Calculated  correction  factors 
are  listed  in  Table  2. 

3.  Simulated  maturation  of  Barney  Creek  sediment  by  MSSVpy 

3.1.  Overall  product  distribution 

A  representative  chromatogram  from  MSSVpy  of  the  BCF  sam¬ 
ple  (300  °C)  is  shown  in  Fig.  1.  The  most  abundant  products  are 
listed  in  Table  3.  Aliphatic  products  from  MSSVpy  were  predomi¬ 
nantly  normal  and  methylated  alkanes  (marked  with  closed  and 
open  circles,  respectively).  Low  molecular  weight  alkanes  (C4-C9) 
displayed  a  higher  proportion  of  methylation.  n-Alkanes  >  Cio 


2 


Fig.  1.  Total  ion  chromatogram  from  MSSVpy  of  the  BCF  sample,  pyrolysed  at  300  °C  for  72  h,  with  no  added  metal  sulfides.  Labelled  peaks  are  listed  in  Table  3. 
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Table  2 

Diagnostic  ions  and  mass  spectral  correction  factors  for  major  PAH  products. 


Compound 

Diagnostic  ion 

Correction  factor 

Phenanthrene 

178 

0.45 

Methylphenanthrenes 

192 

0.34 

Pyrene 

202 

0.37 

Chrysene  +  triphenylene 

228 

0.37 

Benzo[e]pyrene 

252 

0.30 

Benzo[g/n']perylene 

276 

0.46 

displayed  a  unimodal  distribution  with  a  maximum  at  C13,  a  long 
tail  extending  to  C31  and  a  lower  proportion  of  methylated  isomers. 
Low  MW  alkenes,  cycloalkanes  and  ketones  were  present  in  rela¬ 
tively  low  abundances.  The  major  aromatic  products  were  ben¬ 
zenes  and  naphthalenes  with  1-5  methyl  groups.  Phenanthrenes 
and  high  MW  PAHs  with  4-7  aromatic  rings  were  present  in  much 
lower  abundance  (see  Sections  3.2  and  3.3).  Pyrolysis  at  higher 
temperatures  gave  similar  results,  with  the  most  significant  differ¬ 
ence  being  an  increase  in  the  abundance  of  high  MW  PAHs,  consis¬ 
tent  with  the  increase  in  PAH  abundance  observed  across  a  thermal 
gradient  at  HYC  (Chen  et  al.,  2003)  (see  Fig.  2  and  Section  3.3  for 
discussion). 

Solvent  extracted  bitumens  from  the  BCF  contain  abundant  aryl 
isoprenoids  (Brocks  et  al.,  2005;  Brocks  and  Schaeffer,  2008),  which 
are  formed  from  the  breakdown  of  microbial  carotenoids 
(Summons  and  Powell,  1987;  Grice  et  al.,  1996,  1997;  Koopmans 
et  al.,  1996).  These  were  not  detected  in  the  MSSVpy  experiments. 
Pyrolysis  of  long  chain  alkylbenzenes  has  been  shown  to  produce 
methylbenzene  through  cleavage  of  the  alkyl  C-C  bond  in  the  fl- 
position  to  the  aromatic  ring,  a  mechanism  favoured  by  the  high 
stability  of  the  benzyl  radical  intermediate  (Safarik  and  Strausz, 
1997).  Cleavage  of  the  alkyl  chain  from  the  aryl  isoprenoids  pres¬ 
ent  in  the  BCF  bitumen  would  produce  tetramethylbenzenes  with 
a  predominant  1,2, 3, 4  substitution  pattern.  Abundant  1, 2,3,4- 
tetramethylbenzene  was  noted  from  flash  pyrolysis  of  sedimentary 
organics  from  the  Duvernay  Fm.  (Western  Canada  Sedimentary 
Basin)  and  was  attributed  to  the  breakdown  of  di-aromatic  carote¬ 
noids  (Hartgers  et  al.,  1994;  Maslen  et  al.,  2009).  In  the  MSSVpy 
products  from  this  study  the  mono-,  di-  and  trimethylbenzenes 
were  more  abundant  than  the  tetramethylbenzenes.  It  is  likely 
they  were  produced  by  demethylation  and  methyl  shift  reactions 
(e.g.  Szczerba  and  Rospondek,  2010;  Le  Metayer  et  al.,  2014)  result¬ 
ing  in  a  larger  suite  of  polymethylated  products. 

3.2.  Methylated  phenanthrene  ratios 

The  average  values  for  three  commonly  used  methylphenanth- 
rene  (MP)  maturity  parameters  (MPI-1,  MPR  and  MPI-3)  from  the 
replicate  MSSVpy  experiments  are  listed  in  Table  4,  with  the  range 
of  replicate  measurements  in  parentheses.  The  reproducibility  was 
in  most  cases  within  5%.  For  the  pure  BCF  sample,  the  maturity 
parameters  increased  with  pyrolysis  temperature,  consistent  with 
their  behaviour  during  natural  maturation  (Radke  et  al.,  1982a). 
The  average  MPI-1  values  correspond  to  equivalent  vitrinite  reflec¬ 
tance  values  of  0.89  at  300  °C  to  1.16  at  360  °C,  using  the  calibra¬ 
tion  of  Boreham  et  al.  (1988). 

The  calculated  MP  ratios  were  compared  to  those  from  HYC  bit¬ 
umens  reported  by  Williford  et  al.  (201 1 ),  which  are  shown  by  the 
shaded  areas  in  Fig.  3A-C.  The  ratios  using  only  MP  isomers  (MPR 
and  MPI-3)  display  similar  behaviour  in  comparison  to  HYC  bitu¬ 
men.  The  BCF  sample  pyrolysed  at  300  °C  gave  values  slightly 
below  the  minimum  reported  from  HYC.  At  330  °C  the  values  lie 
at  the  middle  of  the  HYC  range,  while  at  360  °C  values  significantly 
exceed  the  maximums  reported  for  HYC.  Unlike  MPR  and  MPI-3, 
the  MPI-1  ratio  includes  unmethylated  phenanthrene  as  well  as 


Table  3 

Major  products  from  MSSVpy  of  the  BCF  sample  at  300  °C  for  72  h,  with  no  added 
metal  sulfides. 


Label 

Compound 

Label 

Compound 

1 

co2 

B4 

C4  benzenes 

2 

Propene 

B5 

C5  benzenes 

3 

2-Methylpropane 

N 

Naphthalene 

4 

3-Methylbut-l-ene 

N1 

Methylnaphthalenes 

5 

Methylpent-2-ene 

N2 

C2  naphthalenes 

6 

Methylcyclopentane 

N3 

C3  naphthalenes 

7 

Methylcyclohexane 

N4 

C4  naphthalenes 

• 

n-Alkanes 

N5 

C5  naphthalenes 

O 

Methyl  alkanes 

P 

Phenanthrene 

B 

Benzene 

PI 

Methylphenanthrenes 

B1 

Toluene 

P2 

C2  phenanthrenes 

B2 

C2  benzenes 

P3 

C3  phenanthrenes 

B3 

C3  benzenes 

P4 

C4  phenanthrenes 

MP  isomers.  Whilst  MPI-1  values  similarly  increased  with  higher 
pyrolysis  temperatures,  they  differed  from  the  range  reported  for 
HYC  (Fig.  3A).  MSSVpy  at  300  °C  and  330  °C  gave  MPI-1  values  that 
were  significantly  lower  than  the  minimum  observed  from  HYC, 
and  only  the  value  at  360  °C  was  within  the  range  of  HYC.  The 
MPI-1  values  from  MSSVpy  were  effectively  lower  than  MPR  and 
MPI-3  in  comparison  to  the  reported  values  from  HYC. 

The  total  abundances  of  the  MP  isomers  relative  to  phenan¬ 
threne  (referred  to  in  this  study  as  MP/P)  is  from  2. 7-3.9  for  all 
MSSVpy  experiments,  and  generally  decreased  slightly  from 
300  °C  to  360  °C  pyrolysis  temperature.  For  HYC  the  MP/P  ratio  is 
from  5.0-9.3  (Williford  et  al.,  2011).  An  increased  proportion  of 
phenanthrene  in  the  MSSVpy  experiments  is  responsible  for  a  low¬ 
ering  of  MPI-1  values.  Pyrolysis  can  promote  the  cleavage  of  alkyl 
groups  from  aromatic  molecules  (Safarik  and  Strausz,  1997),  which 
may  account  for  the  low  abundances  of  MPs  observed  from 
MSSPpy  experiments.  Additionally,  oxidation  of  OM  during  pyroly¬ 
sis  may  result  in  a  lowering  of  the  MPI-1  and  MP/P  ratios.  Phenan¬ 
threne  has  a  higher  average  oxidation  state  of  carbon  than  the  MP 
isomers  (-0.71  compared  to  -0.8).  The  values  of  MPI-1  and  MP/P 
are  therefore  sensitive  to  changes  in  redox  conditions  as  well  as 
temperature,  and  oxidation  of  OM  would  be  expected  to  decrease 
the  MP/P  ratio  and  lower  the  value  of  MPI-1  (Szczerba  and 
Rospondek,  2010;  Dick  et  al.,  2013).  A  reduction  in  MPI-1  corre¬ 
sponding  to  a  shift  from  reducing  to  oxidising  conditions  has  been 
observed  in  the  Kupferschiefer  deposit  (Piittmann  et  al.,  1989;  Sun, 
1998).  However,  the  use  of  glass  beads  to  fill  the  majority  of  void 
space  within  the  tube  essentially  renders  the  MSSVpy  technique 
free  of  oxidation  artefacts  (Horsfield  et  al.,  2015).  Indeed  the 
observed  increase  of  MPI-1  with  pyrolysis  temperature  (Fig.  3A) 
shows  that  temperature  is  the  major  control  on  MPI-1  in  these 
experiments. 

Alternatively  an  additional  mechanism  that  was  not  present  in 
the  MSSVpy  experiments  may  have  produced  MPs  at  HYC.  One 
potential  explanation  is  the  migration  of  aromatic  hydrocarbons 
to  HYC  with  the  mineralising  fluid.  Williford  et  al.  (2011)  found 
that  PAHs  from  HYC  have  similar  <513C  to  kerogen  from  the  under¬ 
lying  Wollogorang  Fm.  and  proposed  that  PAHs  were  generated  in 
the  Wollogorang  Fm.  by  hydrothermal  alteration  and  deposited  in 
the  BCF  by  the  mineralising  fluid.  It  is  unclear  however  if  these 
PAHs  would  have  a  greater  ratio  of  MP/P  than  those  produced  by 
MSSVpy  in  this  study.  The  Wollogorang  Fm.  has  experienced 
extensive  metamorphism  from  volcanic  intrusions  (Donnelly  and 
Jackson,  1988),  which  has  likely  resulted  in  the  destruction  of  all 
indigenous  organic  matter  (Brocks  et  al.,  2008).  The  composition 
of  indigenous  OM  from  this  formation  therefore  remains  uncertain 
and  the  likelihood  of  this  explanation  cannot  presently  be 
evaluated. 
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Fig.  2.  Summed  m/z  202  +  228  +  252  +  276  chromatograms  from  MSSVpy  of  the  BCF 
sample  at  each  pyrolysis  temperature,  with  no  added  metal  sulfides.  Pyrolysis 
temperature  are  (A)  300  °C,  (B)  330  °C,  (C)  360  °C.  Fla:  fluoranthene,  Pyr:  pyrene, 
BaAn:  benzo[a]anthracene,  Chry:  chrysene,  Tpn:  triphenylene,  BePy:  benzo[e]pyr- 
ene,  BaPy:  benzo[a]pyrene,  BPery:  benzo[gfii]perylene. 


3.3.  High  MW  PAHs 

PAHs  with  4-6  aromatic  rings  (hereby  collectively  referred  to  as 
high  MW  PAHs)  were  detected  at  low  abundance  in  all  MSSVpy 
experiments.  Summed  m/z  202  +  228  +  252  +  276  chromatograms 
for  the  three  pyrolysis  temperatures  are  shown  in  Fig.  2.  At 
300  °C  the  four-ring  PAHs  fluoranthene  and  pyrene  have  lower 
abundance  than  other  peaks  in  the  m/z  202  chromatogram  (mainly 
dimethyl-  and  trimethylphenanthrenes),  while  the  five  and  six- 
ring  PAHs  are  present  only  in  trace  amounts.  At  330  °C  and 
360  °C  the  PAHs  become  far  more  prominent,  implying  increased 
production  of  high  MW  PAHs  at  higher  pyrolysis  temperature.  This 


Table  4 

PAH  ratios  measured  from  MSSVpy  analyses.  Values  are  reported  as  the  average  of 
duplicate  measurements,  with  the  range  of  measurement  in  parentheses. 


300  °C 

330  °C 

360  °C 

BCF 

MPI-1 

0.96  (0.02) 

1.07  (0.05) 

1.34(0.01) 

MPR 

1.07  (0.06) 

1.24  (0.04) 

1.69  (0.06) 

MPI-3 

1.05  (0.03) 

1.15  (0.07) 

1.58  (0.01) 

MP/P 

3.21  (0.01) 

3.51  (0.04) 

3.39  (0.02) 

PAH/P 

0.45  (0.01) 

0.77  (0.04) 

0.91  (0.02) 

BCF  +  5%  PbS 

MPI-1 

1.04  (0.01) 

0.98  (0.08) 

1.17  (0.03) 

MPR 

1.15  (0.03) 

1.23  (0.03) 

1.52  (0.04) 

MPI-3 

1.11  (0.04) 

1.08  (0.06) 

1.34  (0.01) 

MP/P 

3.47  (0.13) 

3.24  (1.07) 

3.24(0.16) 

PAH/P 

0.37  (0.05) 

0.57  (0.29) 

0.51  (0.31) 

BCF  +  20%  PbS 

MPI-1 

1.14  (0.06) 

1.02  (0.04) 

1.15  (0.01) 

MPR 

1.25  (0.05) 

1.26  (0.03) 

1.42  (0.04) 

MPI-3 

1.21  (0.03) 

1.09  (0.08) 

1.28  (0.01) 

MP/P 

3.75  (0.31) 

3.47  (0.18) 

3.42  (0.05) 

PAH/P 

0.24  (0.29) 

0.43  (0.03) 

0.69  (0.03) 

BCF  +  5%  ZnS 

MPI-1 

1.06  (0.12) 

1.06  (0.01) 

1.16  (0.01) 

MPR 

1.21  (0.11) 

1.24  (0.05) 

1.49  (0.13) 

MPI-3 

1.20(0.03) 

1.13  (0.01) 

1.38  (0.11) 

MP/P 

3.22  (0.72) 

3.55  (0.09) 

3.08  (0.42) 

PAH/P 

0.29  (0.24) 

0.70  (0.16) 

0.53  (0.19) 

BCF +  20%  ZnS 

MPI-1 

1.03  (0.12) 

1.03  (0.01) 

1.15  (0.01) 

MPR 

1.12  (0.03) 

1.20  (0.01) 

1.43  (0.01) 

MPI-3 

1.16  (0.01) 

1.17  (0.01) 

1.36  (0.01) 

MP/P 

3.19  (0.92) 

3.09  (0.09) 

3.07  (0.01) 

PAH/P 

0.30  (0.3) 

0.35  (0.08) 

0.57  (0.01) 

BCF  +  5%  Py 

MPI-1 

1.21  (0.08) 

1.14  (0.02) 

1.30  (0.01) 

MPR 

1.31  (0.03) 

1.22  (0.01) 

1.51  (0.01) 

MPI-3 

1.29  (0.02) 

1.21  (0.04) 

1.47  (0.01) 

MP/P 

3.85  (0.85) 

3.78  (0.13) 

3.57  (0.01) 

PAH/P 

0.34  (0.19) 

0.78  (0.13) 

0.57  (0.01) 

BCF  +  20%  Py 

MPI-1 

1.17  (0.04) 

1.02  (0.08) 

1.26  (0.01) 

MPR 

1.30(0.08) 

1.15  (0.01) 

1.52  (0.01) 

MPI-3 

1.31  (0.15) 

1.13  (0.08) 

1.49  (0.01) 

MP/P 

3.45  (1.01) 

3.17  (0.2) 

3.21  (0.01) 

PAH/P 

0.36  (0.36) 

0.73  (0.2) 

0.54  (0.01) 

is  consistent  with  observations  across  a  thermal  gradient  at  HYC 
(Chen  et  al.,  2003;  Williford  et  al.,  2011). 

The  summed  abundance  of  the  high  MW  PAHs  fluoranthene, 
pyrene,  benzo[a]anthracene,  chrysene,  triphenylene,  benzo[e] 
pyrene,  benzo[a]pyrene  and  benzo[g/ii]perylene  relative  to  the 
abundance  of  phenanthrene  is  defined  here  as  PAH/P.  This  ratio 
increases  with  pyrolysis  temperature,  approximately  doubling 
from  300  °C  to  360  °C  (see  Table  4  and  Fig.  3D).  Phenanthrene 
was  reported  to  be  totally  stable  during  hydrous  pyrolysis  at 
330  °C  for  up  to  17  days  (McCollom  et  al.,  1999)  and  even  persisted 
after  pyrolysis  at  400  °C  for  30  days  (Price,  1993),  therefore 
destruction  of  the  phenanthrene  ring  system  is  not  expected  to 
occur  at  the  conditions  used  in  this  study.  The  increased  PAH/P 
ratio  is  attributed  to  the  increased  production  of  high  MW  PAHs 
at  higher  pyrolysis  temperatures. 

Bitumen  isolated  from  the  HYC  deposit  contains  a  high  propor¬ 
tion  of  high  MW  PAHs,  with  the  samples  from  Williford  et  al. 
(2011)  displaying  PAH/P  ratios  from  2. 2-7.9.  In  contrast  all 
MSSVpy  experiments  showed  PAH/P  ratios  of<  1.0  (Fig.  3D)  due 
to  the  greatly  reduced  abundances  of  high  MW  PAHs.  It  has  been 
proposed  that  high  MW  PAHs  in  hydrothermal  systems  are  formed 
by  the  condensation  of  low  MW  organic  precursors  (Kawka  and 
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■  BCF  A  5%  PbS  ▲  20%  PbS  0  5%ZnS 
•  20%  ZnS  O  5%  Py  ♦20%Py 


Fig.  3.  Selected  PAH  ratios  from  MSSVpy.  (A)  methylphenanthrene  index  1  (MPI-1),  (B)  methylphenanthrene  ratio  (MPR),  (C)  methylphenanthrene  index  3  (MPI-3)  and  (4) 
ratio  of  four  to  six  ring  PAHs  over  phenanthrene  (PAH/P).  Data  points  are  the  average  of  duplicate  measurements,  error  bars  show  the  range  of  measurement.  Shaded  areas  in 
A  to  C  represent  the  range  from  the  HYC  deposit  reported  by  Williford  et  al.  (201 1 ).  The  range  of  PAH/P  reported  by  Williford  et  al.  (201 1 )  is  2.2-7. 9. 


Simoneit,  1990;  Fetzer  et  al.,  1996)  through  the  addition  of  C2  or 
C4H2  units  to  existing  aromatic  structures  (Stein,  1978).  There 
would  likely  be  sufficient  organic  precursors  present  in  the  BCF 
sediment  to  form  significant  amounts  of  high  MW  PAHs,  suggest¬ 
ing  that  the  low  abundances  of  these  compounds  were  due  to 
the  thermal  conditions  used  during  MSSVpy  treatments.  Higher 
temperatures  or  longer  term  experiments  may  be  required  to  gen¬ 
erate  large  abundances  of  high  MW  PAHs.  Chen  et  al.  (2003)  pro¬ 
posed  that  high  MW  PAHs  were  concentrated  at  HYC  by 
precipitation  from  hydrothermal  fluid,  due  to  the  lower  aqueous 
solubility  of  heavy  aromatics  compared  to  two-  and  three-ringed 
PAHs  (Kawka  and  Simoneit,  1990).  Such  a  mechanism  may  contrib¬ 
ute  to  the  increased  abundances  of  high  MW  PAHs  at  HYC  com¬ 
pared  to  the  MSSVpy  experiments. 

4.  Effects  of  metal  sulfides  on  simulated  thermal  maturation 

4.3.  Methylated  phenanthrene  ratios 

The  addition  of  PbS,  ZnS  and  pyrite  had  a  marked  effect  on  all 
MP  ratios  (Fig.  3A-C).  At  300  °C  pyrolysis  temperature  all  ratios 
were  higher  than  the  pure  BCF  sample.  At  330  °C  the  ratios  were 
roughly  equal  with  and  without  added  metal  sulfides,  while  at 
360  °C  the  ratios  were  notably  lower  with  added  metal  sulfides. 
Some  samples  also  displayed  a  slight  decrease  in  MP  ratios  from 
300  °C  to  330  °C,  suggesting  an  inversion  of  the  maturity  parame¬ 
ters  in  this  temperature  range.  MP  maturity  ratios  from  experi¬ 
ments  with  each  sulfide  mineral,  and  with  5  wt%  and  20  wt% 
sulfide  addition,  were  all  generally  within  uncertainty. 

MP  maturity  ratios  may  be  influenced  by  the  preferential 
adsorption  of  certain  isomers  onto  the  sulfide  mineral  surfaces. 
In  a  previous  study  of  the  HYC  deposit  it  was  found  that  the 
p-MP  isomers  are  more  readily  extracted  with  organic  solvent  than 
the  a  isomers  (Holman  et  al.,  2012).  Preferential  retention  of  a-MP 


isomers  by  mineral  surfaces  would  result  in  a  decrease  in  apparent 
maturity.  It  is  likely  however  that  adsorption  effects  played  only  a 
minor  role  in  the  observed  behaviour  of  the  MP  ratios  in  this  study. 
Strong  organic  adsorption  is  generally  observed  on  clay  minerals 
(Lahav  and  Chang,  1976;  Cleaves  et  al.,  2012),  not  metal  sulfides. 
Adsorption  effects  are  expected  to  be  proportional  to  mineral 
surface  area  (Cleaves  et  al.,  2012),  so  it  would  be  expected  that 
experiments  with  20  wt%  sulfide  minerals  added  would  display 
stronger  effects  than  those  with  5  wt%  sulfides.  Fig.  3  shows  this 
is  not  the  case,  with  little  difference  being  observed  between 
5  wt%  and  20  wt%  sulfides  in  most  cases.  This  argues  against  the 
strong  influence  of  adsorption  effects. 

Alternatively  the  metal  sulfides  may  be  acting  to  catalyse  or 
retard  the  isomerisation  of  MPs.  Transition  metal  sulfides  are 
known  to  act  as  surface  catalysts  for  a  number  of  sedimentary 
organic  reactions  (Wachtershauser,  1988;  Cody  et  al.,  2004).  Dur¬ 
ing  thermal  maturation  MPs  undergo  a  range  of  methyl-shift  and 
transmethylation  reactions  (Szczerba  and  Rospondek,  2010)  which 
may  be  retarded  by  the  presence  of  metal  sulfides.  The  overall 
effect  on  the  MP  maturity  parameters  will  be  determined  by  the 
kinetics  of  each  individual  isomerisation  reaction  and  will  require 
more  detailed  studies  to  fully  elucidate. 

4.2.  High  MW  PAHs 

The  addition  of  sulfide  minerals  caused  a  decrease  in  the  pro¬ 
duction  of  high  MW  PAHs  (Fig.  3D).  This  is  especially  apparent  at 
higher  pyrolysis  temperatures,  where  the  average  PAH/P  ratio 
decreased  from  0.91  for  the  pure  BCF  sample  to  between  0.51 
and  0.69  when  sulfide  minerals  were  added.  Further  interpretation 
of  this  relationship  is  limited  by  the  high  degree  of  uncertainty  in 
the  analysis  of  several  samples,  with  the  PAH/P  ratio  of  some  rep¬ 
licates  varying  by  >  50%.  Fluoranthene  and  chrysene/triphenylene 
are  both  significantly  affected  by  co-eluting  peaks  even  when 
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integrated  using  the  m/z  202  and  228  chromatograms,  while  five 
and  six  ring  PAHs  were  of  such  low  abundance,  especially  at 
300  °C  pyrolysis  temperature  (Fig.  2A),  that  they  were  difficult  to 
integrate  above  the  background  noise.  Nevertheless  the  ratio  of 
PAH/P  for  all  MSSVpy  experiments  with  added  sulfide  minerals 
was  lower  than  those  for  the  pure  BCF  sample.  The  reason  for  this 
reduction  is  currently  unclear.  The  metal  sulfide  surfaces  may  inhi¬ 
bit  the  addition  reactions  by  which  high  MW  PAHs  are  assembled 
from  smaller  aromatic  units.  As  discussed  in  Section  3.2,  the  MPI-1 
and  PAH/P  ratios  are  sensitive  to  changes  in  redox  conditions.  Sul¬ 
fide  minerals  are  redox  active  and  thus  may  influence  the  redox 
conditions  of  the  MSSVpy  experiments,  however  oxidation  during 
pyrolysis  is  negligible  due  to  the  very  low  void  space  within  the 
MSSV  tubes  (Horsfield  et  al.,  2015). 

5.  Conclusions 

MSSVpy  has  been  successfully  used  to  simulate  thermal  matura¬ 
tion  of  OM  from  the  BCF.  MSSVpy  of  thermally  immature  sediment 
produced  a  distribution  of  hydrocarbon  products  including  n- 
alkanes,  two  and  three-ring  PAHs  and  their  methylated  isomers, 
typical  of  thermally  mature  organic  sediments  associated  with  min¬ 
eralised  regions  of  HYC.  Methylphenanthrene  maturity  parameters 
measured  from  MSSVpy  of  BCF  OM  fell  within  the  range  previously 
reported  from  HYC,  although  the  MPI-1  ratio  was  comparatively 
lower  due  to  a  greater  proportion  of  unmethylated  phenanthrene. 
PAHs  with  4-6  aromatic  rings  were  produced  in  increasing  amounts 
at  higher  pyrolysis  temperatures,  but  their  relative  abundance  does 
not  approach  the  levels  reported  at  HYC.  The  temperature  and  time 
conditions  used  in  this  study  may  have  been  insufficient  to  form 
these  compounds  in  high  abundances.  The  increased  maturity 
parameters  and  the  production  of  high  MW  PAHs  were  similar  to 
those  attributed  to  the  hydrothermal  alteration  of  sedimentary 
OM.  The  addition  of  metal  sulfides  consistently  affected  the  values 
of  MP  maturity  ratios,  with  an  increase  at  300  °C  pyrolysis  temper¬ 
ature  and  a  decrease  at  360  °C.  It  is  hypothesised  that  MP  isomerisa¬ 
tion  reactions  may  be  moderated  by  the  sulfide  minerals. 
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